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Abstract 

In planar cellular systems m n denotes the average sidedness of a 
cell neighboring an n-sided cell. Aboav's empirical law states that nm n 
is linear in n. A downward curvature is nevertheless observed in the 
numerical nm n data of the Random Voronoi Froth. The exact large-n 
expansion of m n obtained in the present work, viz. m n = 4+3(7r/n) 2 + 
. . ., explains this curvature. Its inverse square root dependence on n 
sets a new theoretical paradigm. Similar curved behavior may be 
expected, and must indeed be looked for, in experimental data of 
sufficiently high resolution. We argue that it occurs, in particular, in 
diffusion-limited colloidal aggregation on the basis of recent simulation 
data due to Fernandez- Toledano et al. [Phys. Rev. E 71, 041401 
(2005)] and earlier experimental results by Earnshaw and Robinson 
[Phys. Rev. Lett. 72, 3682 (1994)]. 
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1 Introduction 



1.1 General 

In nature, planar cellular systems come in a wide variety. They include 
biological tissues [HEIIS], polycrystals cells formed by particles trapped 
at a water/air interface pj, cells in surface-tension driven Benard convection 
[0], in two-dimensional soap froths i7j, and in magnetic liquid froths [8 . In 
other systems cellular structure may appear when the data are subjected to 
the Voronoi construction [5J. Examples of these are hard disks on an air 
table [Hi], a binary liquid during late stage coarsening [TTJ, two-dimensional 
colloidal aggregation O EH E], nanostructured cellular layers [TH], and 
studies of two-dimensional melting [16 . 

Two empirical rules play a prominent role in studies of planar cellular 
systems: Lewis' law and Aboav's law. Both are statements about the statis- 
tics of a cell's most conspicuous properties, viz. its area and its number of 
sides. Lewis' law [TJ says that the average area A n of an n-sided cell increases 
with n as 

A n = —(n-n ), (1.1) 

where clq and n are constants and A is the areal cell density. In the present 
work we are interested in the second one of these laws, formulated by Aboav 
jl], who noticed that many-sided cells tend to have few-sided neighbors and 
vice versa. He expressed this correlation in terms of the average m n of the 
number of sides of a cell that neighbors an n-sided cell. Aboav's law, also 
called the Aboav- Weaire [H] law, asserts that 

b . . 

m n = 6-a-\ — , (1.2) 
n 

where a and b are numerical constants. This law, which expresses an intu- 
itively plausible trend with n, is in widespread use El El El El El EDI Ell 
IT2*1 [T31 EH] in the analysis of experimental data on cellular structures. One 
usually plots nm n = (6 — a)n + b as a function of n and often refers to either 
this relationship or equation (jl.2j) as the linear law. 

The surprising fact is that for many of these cellular systems, in spite of all 
their diversity, Aboav's linear law appears to hold with good accuracy in most 
of the experimentally accessible range, which runs from n = 3 to n typically 
between 9 and 12. The experimental values of the Aboav parameters, listed 
e.g. in reference [TS], are typically in the range 0.7 < a < 1.5 and 5.7 < 
b < 8.5 [121 . Attempts to explain the linearity with n and the numerical 
values of a and b from first principles, i.e., on the basis of a microscopic 
geometrical model, have been unsuccessful. The various derivations of ()1.2j) 
that one does find in the literature all involve approximations (usually of the 
mean-field kind; see e.g. reference |20j) or hypotheses whose general validity 
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is subject to caution (we consider the 'maximum entropy method' [21] to be 
in this class). The dominant view today is probably that ()1.2|) is not exact, 
but merely a good approximation to some unknown 'true' curve, which need 
not be the same for all cellular systems. 

It has been realized, in particular by Le Caer and Delannay [22], that 
knowledge about the large-n behavior ofm n will constrain the law for m n , for 
example by establishing the regime of validity of its linearity or by putting 
limits on the numerical values of a and b. The question of the large-n behavior 
of m n is however a difficult one and has received little attention since. The 
'derivations' referred to above generally lead to (jl.2j) without providing any 
restrictions on its range of validity. In the present work we return to the 
large-n behavior. We analyze it within the context of the Random Voronoi 
Froth and consider the implications of our conclusion for the interpretation 
of simulations and experiments. 

1.2 The Random Voronoi Froth (RVF) 

We consider the 'Random Voronoi Froth' (RVF), in the mathematical lit- 
erature rather called the 'Poisson- Voronoi tessellation'. It is obtained by 
constructing the Voronoi cells [9 of a planar set of randomly and uniformly 
distributed point centers, for convenience often called 'seeds' (but without 
the implication that they are material). The statistical properties of the 
RVF have been studied more than those of any other microscopic model of 
cellular structure. Reference P compiles a large body of analytic results and 
numerical tables. For the role and place of the RVF amidst other models of 
cellular systems one may consult, e.g., Rivier [22] or Schliecker |24j . 

Analytically, a first-principle derivation of an exact expression for m n in 
the RVF (or, for that matter, in any other geometrical model) is still lack- 
ing. The reason is that the calculation of this quantity is in the notoriously 
difficult class of many-body problems. Numerically, however, the m n values 
of the RVF are known with considerable precision fI5\ l2"o1 12~71 128j . The best 
data come from Monte Carlo simulations by Brakke [2H] , whose results for m n 
have a four digit accuracy in the range 4 < n < 9. In figure show these 
data as well as those by Boots and Murdoch j2H], represented in the usual 
way in an nm n versus n plot. The dashed straight line is Aboav's law (|1.2|) 
with a = 0.75 and b = 5.76. Although this law provides what appears as a 
very good approximation, the data points nevertheless exhibit an extremely 
small but distinctive downward curvature and the narrowness of their error 
bars rules out Aboav's linear fit. We have computed the second derivative 
(nm n )" = (n — l)m n -\ — 2nm n + (n + l)m n+ i from Brakke's [2H] numerical 
data and plotted it in figure H The ratio —(nm n )"/(nm n ), which may be 
taken as a measure of the curvature, is of the order of only 1/250 = 0.004. 
Hence a successful theory of the RVF should explain the existence as well as 
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the sign and the smallness of the curvature. So far, in order to accommodate 
this departure from linearity, only ad hoc alternatives to (jl.2J) have been 
proposed (see, e.g., the discussions in j2H] and in section 5.3.3 of |9J. 

We pursue here the approach initiated in reference |2H|, where the meth- 
ods of statistical mechanics were brought to bear on planar Voronoi tessel- 
lations. That work has opened the possibility of an expansion in powers of 
for all quantities of interest related to an n-sided Voronoi cell. This 
initially led to the asymptotic large-n expansion of the probability p n for an 
arbitrarily chosen cell to have n sides. An immediate further result was the 
proof that the RVF obeys Lewis' law (jl.ljl asymptotically for n — > oo with 
coefficient a$ = j. But whereas Lewis' law refers to a single cell, Aboav's 
deals with the intrinsically more difficult problem of cell-cell correlations; it 
therefore requires the separate study that we present here. 

Our key result, equation (|2.2|) of section 12.11 is an exact expression for 
the asymptotic behavior of m n as n becomes large. The inverse square root 
decay appearing in that expression explains the small deviations from linear- 
ity that are observed. As a corollary we obtain in section |2~2*1 the asymptotic 
expression for the total number, to be called Kn , of second-neighbor cells 
surrounding a central n-sided cell. In section |3] we discuss these results. We 
compare them to the numerical RVF data, as well as to data on diffusion 
limited colloidal aggregation coming from both experiments and simulations. 
Whereas our asymptotic expansion constitutes an exact result based on first 
principles, in section 0] we take a pragmatic attitude and try to find the best 
two- and three-parameter fits for m n that incorporate this large n behavior. 
In section |5] we speculate briefly on modifications that our theory may un- 
dergo in the case of cellular systems not of the RVF type. In section El we 
conclude. Finally, in appendix A we prove a theorem concerning the statis- 
tics of randomly distributed points in a half-plane that is indispensible for 
the derivation of the main result. 

2 Two-cell correlations 
2.1 First-neighbor correlation 

We will determine analytically the large-n behavior of m n in a RVF. Two 
steps are required to arrive at this result. The first step of the derivation 
is based on the geometrical properties of the large n-sided cell; these have 
all been determined quantitatively in earlier work where they required 
considerable mathematical effort. The arguments of the first step will be 
presented with the aid of a figure, but we stress that they are nonetheless 
nontrivial and dictated by compelling logic; they constitute the conceptually 
new part of this investigation. The second step is a problem in statistics 
whose solution, although somewhat lengthy, can be obtained by standard 
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Figure 1: The cells neighboring an n-sided cell have together a total average 
of nm n sides. Open circles: simulation values of nm n according to Boots and 
Murdoch j2^j. Dots: simulation values of nm n according to Brakke [2Sj. Dashed 
line: Aboav's linear law nm n = (6 — a)n + b with a = 0.75 and b = 5.76. Open 
squares: the second derivative (nm n )" = (n — l)m n _i — 2nm n + (n + l)m n+ i, 
constructed from Brakke's [2S] i^n data and multiplied by —250 in order to be 
visible on the scale of the figure. Error bars are shown for the dot and open square 
data; where absent, they are smaller than the data point symbols. 



methods. We present it in appendix 1X1 

We begin by considering an n-sided cell with n very large. Cells with very 
many sides are extremely rare, but if one occurs, then its environment must 
look as depicted in figure El where the n-sided cell of a "central" seed at O is 
surrounded by n strongly elongated first neighbor cells containing seeds Fj. 
Independent evidence for such a geometry comes from work by Lauritsen et 
al. [30J, who Monte Carlo simulated a Hamiltonian favoring the appearance 
of many-sided cells. In figure El four different length scales [29J play a role, 
each proportional to its own characteristic power of n: 

(i) The perimeter T of the central n-sided cell typically runs within an 
annulus of center O, of radius R c = (n/47rA)2, and of width of order 1. 
Hence, provided it remains smooth, the perimeter tends for n — > oo towards 
a circle of radius R c . Similarly, the first neighbor seeds Fi will be on a circle 
of radius 2R C . 

(ii) Sufficient smoothness of the perimeter Tq in the limit n — > oo is 
guaranteed by the following property: locally the radial coordinates of the 
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Figure 2: Schematic picture of the environment of an n-sided Voronoi cell with 
n ~ 100 around a seed at O. The first and second neighbor cells have their seeds 
at F\, F2, -F3, . . . and at Si, S2, S3, . . ., respectively. All solid line segments separate 
Voronoi cells. Among these, the heavy solid line Tq is the perimeter of the n-sided 
cell, which is close to circular. The heavy solid line Ti separates the first from the 
second neighbors. Both Tq and T± are piecewise linear on a scale n~2. On the 
scale of order 1 the incipient piecewise parabolic structure of Ti is discernible. 
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individual vertices of Tq have rms deviations of order n~2 with respect to the 
locally averaged radial coordinate. A similar statement holds for the curve, 
not drawn in the figure that links the successive first neighbors. 

(iii) The n vertices on the perimeter of the central cell have a line density 
Pvort = n/(2nR c ) = (nX/ir)?. Consequently, the first neighbor seeds Fj have 
a line density |p VC rt- 

(iv) In the region to the right of the heavy solid line r 1; which is occupied 
by second and further neighbors, the seed density keeps its "background" 
value A, i.e. is of order n° [HTj . 

On this picture we base the following line of arguments. For n — ► oo 
the vertex line density p ver t = {nX/n)^ tends to infinity. In spite of this the 
areal density A of the seeds to the right of Ti stays of order n°. It follows 
that the central cell will have ~ second neighbor cells. Since there are n 
first neighbors, each second neighbor cell Sj must be adjacent to ~ first 
neighbor cells Fi (where we call the cells by the names of their seeds and 
where the symbol ~ denotes asymptotic proportionality). Figure El shows 
that under these geometrical constraints each first neighbor is most likely 
to have itself four neighbors, viz. the central n-sided cell, a single second 
neighbor cell, and two other first neighbors, Fi-i and i^+j. 

We now focus on the exceptional Fi that have five neighbors due to their 
being adjacent to two second neighbors Sj and Sj+i. An example is the cell 
marked F\ in figure El which is adjacent to both S\ and S 2 - Let us denote 
by /g the fraction of first neighbors that are five-sided. In view of the scaling 
relations that precede we expect that /g = cn~2 + . . . , where c is a numerical 
coefficient and the dots indicate terms of higher order in n~2. Any six- and 
higher-sided Fi will contribute only to these dot terms. Hence we have 

m n = 4 (l-/ 5 ) + 5/ 5 

= 4 + cn"5 + . . . . (2.1) 

It turns out to be possible to determine the constant c. To that end we 
consider the solid line T\ in figure which separates the central seed's first 
from its second neighbors. In the large-n limit the curve linking the Fi 
becomes a circle which may locally be replaced with a straight line. In 
figure H3 this straight line is represented by the x axis and the region of space 
containing the second and further neighbors by the half-plane y > 0. Seeds 
are uniformly distributed in the upper half-plane with the background density 
A. Since the first neighbors Fi are dense on the x axis, the curve Ti divides 
the half-plane y > into a lower part of points closer to the x axis than 
to any of the seeds, and its complement. Hence the function y = Ti(x) is 
piecewise parabolic; its incipient parabolic segments are discernable in figure 
121 To each cusp of Ti(x) corresponds a five-sided first neighbor cell. Let 
Pcusp t> e the density on the x axis of the abscissae of the cusps of T^x). Then 
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Figure 3: The x-axis represents a continuum of first neighbors. The upper half- 
plane is randomly filled with seeds of uniform density A (here A = 1). The region 
of the half-plane which is closer to the x axis than to any of the seeds is separated 
from its complement by the piecewise parabolic curve y = ri(x). The abscissae of 
the cusps of Ti are shown in appendix A to have a density of §A5 on the x axis. 
Full and open squares represent second and further neighbors, respectively. Note 
that the horizontal and vertical scales are different. 
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it is clear that /s = 2p cnsp / p vcrt . The determination of p CU s P for given seed 
density A in the upper half plane is a well-defined problem in statistics. A 
somewhat lengthy calculation yields p cnsp = |Aa (see appendix E} , in which 
the coefficient | is essential. Using the expression for p vcrt found above we 
therefore have that /s = 2(ir/n\)^ x 3X^/2 = 3(ir/n)2, whence c = 37T2. 
Substituting this in ()2.1|) we conclude that for the Random Voronoi Froth 
m n is given by 

m n = 4 + 3 A /- + ..., (n->oo), (2.2) 
V n 

where the dots stand for higher order terms in n~2. This exact asymptotic 
formula for the first-neighbor correlation of Poisson- Voronoi cells is the key 
result of this work. It represents a new paradigm for the behavior of the 
cell-cell correlation. We defer all further comments to section El 



2.2 Second-neighbor correlation 

It is easy to obtain a corollary involving the correlation between a central 
cell and its second neighbor cells in the RVF. In a wider context concentric 
layers of cells around a central one have been the subject of various investiga- 
tions in recent years [32] ■ We denote by the average number of second 
neighbors of an n-sided cell. The preceding result allows us to determine 

(2,) 

almost immediately the expression for KX\ again in the limit of asymptot- 
ically large n. We appeal once more to figure El For n — > oo the circle of 
seeds Fi has a circumference 4irR c = 2(mr/X)2, The total number N cusp of 
cusps of Ti along this circle is equal to p CU s P times the circumference, whence 
N cusp = 3 (7m) 2. Some caution must be exercised at this point. The number 
of second-neighbor cells can at most be equal to N cusp ; however, it will be 
slightly smaller. The reason is that the common border between a given sec- 
ond neighbor and the set of first neighbors may consist of more than a single 
parabolic segment. We therefore have K n = r/N cusp where the reduction 
factor 1] < 1 is a geometrically defined mathematical constant. The result is 
that 

= 377(7m)2 + ..., (rwoo). (2.3) 

We have not seen a way to find a simple analytic expression for rj; from 
a simulation in which we generated more than 2000 cusps we obtained the 
estimate 

77 = 0.95 ±0.01. (2.4) 

Equation ()2.3j) together with (|2.4jl constitutes the exact asymptotic formula 
for the total number of second neighbors of a central n-sided cell in a Random 
Voronoi Froth. 
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3 Discussion 

We now discuss the significance of equation (j2.2jl . This equation is first of all 
a statement of principle: in the RVF (and hence in all theoretical models or 
experimental systems systems for which it is relevant), as the accuracy and 
the range of the data increase, the two-cell correlation m n should be seen 
to follow an n~2 law and asymptotically approach = 4 [331 - Equation 
(12. 2|) therefore rules out the possibility of the existence of a general proof of 
Aboav's law for an unrestricted range of n. Aboav's law keeps its meaning, 
however, as a good linear approximation to the available data in the window 
accessible to simulations or experiments. 

3.1 Comparison to RVF simulation data 

Earnshaw and Robinson ^2j, following a remark for which they credit Weaire, 
have emphasized that in order to detect small deviations from Aboav's law 
it is essential to plot the m n data as a function of n _1 . We fully concur with 
them and have plotted in this way in figure HJthe same data that were shown 
in figure H in the new variable n^ 1 Aboav's law remains linear and the 
negative curvature of the simulation data now appears clearly. In addition 
we have represented by a solid line in figure 0] the first two terms, m n = 

4 + 3(7r/n)a, of equation (|2.2j) . We discuss separately two aspects of the 
theory, viz. its prediction of the curvature effect and its numerical accuracy. 

Curvature. - We confront our calculation of the asymptotic behavior of 
m n with the RVF simulation data. It is visually apparent from figure H] that 
the simulation data have a curvature very similar to that of the theoretical 
curve. Quantitatively, we obtain from (J2.2)) analytically 

(nm n )" 37r^ 
nm n 4 n f [4 + 3(71-/71)*] ' 

which ranges from 0.006 to 0.001 in the interval 4 < n < 8 and hence is of the 
same order of magnitude as the value 0.004 found from figure^ We conclude 
that the asymptotic result (|2.2|) not only explains why the m n data should 
be curved, but correctly predicts the sign and magnitude of the curvature. 

Numerical accuracy. - Equation (|2.2|) is the beginning of an asymptotic 
series in powers of rT^. Although they may, such series are not guaranteed 
to produce accurate numerical values when the expansion variable is finite. 
Figures 121 (which is for n « 100) and El (for n = 00), used in our derivation, 
refer to seed configurations that are very far away indeed from those observed 
in either simulations or experiments. Therefore, one might be uninclined 
a priori to expect of this expansion a high degree of numerical accuracy. 
However, figure 0] shows the contrary. When compared to the numerical 
RVF values, our asymptotic result (|2.2j) . limited to its first two terms, is 
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Figure 4: Average sidedness m n of a neighbor of an n-sided cell plotted against 
n _1 . Open circles and closed dots: same data as in figure ^ but in this plot the 
curvature appears more clearly. Dashed straight line: Aboav's law m n = a + bn^ 1 
with a and b as in figure ^ Solid line: first two terms, m n = 4 + 3(7r/n) 5 , of the 
exact asymptotic series (|2.2|l . Error bars are shown for the dotted data; where 
absent, they are smaller than the data points. 

quite good: it is only slightly more off than Aboav's linear fit. This gives 
confidence in the expansion not only as an asymptotic constraint on the m n 
curve, but also as a practical tool for estimating this and other correlations. 
We will present further numerical considerations in section HJ 

3.2 Comparison to diffusion limited colloidal aggrega- 
tion (DLCA) 

The RVF is an idealized model which appears in the discussion and analy- 
sis of various naturally occurring cellular systems We now investigate 
the relevance of the RVF cell-cell correlations calculated in section El for 
a real physical system, viz. diffusion limited colloidal aggregation (DLCA) 
in two-dimensional suspensions. In the experiments one monitors the slow 
aggregation of particles trapped at the air/water interface and initially ran- 
domly distributed. In the early stage of this process isolated clusters appear 
and video images of the system taken after some time t are analyzed in terms 
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of the Voronoi cells constructed around the centers of mass of these clusters. 
Simulation methods developed for such aggregating systems are known [T3] 
to agree well the experimental studies and we will discuss the simulations 
first. 

Recently detailed Brownian dynamics simulations were performed by 
Fernandez- Toledano et al. [34J. Their m n versus 1/n data have been rep- 
resented in our figure The authors fit the m n data by Aboav's linear 
law (upper thin straight line in figure EJ), but their data points depart from 
linearity for n > 9 in an even more pronounced way than those of figure 
13] Fernandez- Toledano et al. attribute these deviations to statistical uncer- 
tainties. Instead, we believe that their simulations were accurate enough for 
them to actually see the curvature predicted by our asymptotic RVF theory, 
represented in figure El by the solid curve. It predicts a curvature of the same 
order of magnitude as observed in the data. While both Aboav's law and our 
curve are compatible with the error bars, our theory provides an explanation 
whereas the linear fit does not. The authors of reference [33] fit their Krp 
data, shown in figure El with a linear curve based on Aboav's law plus some 
additional hypotheses. The data show, however, a slight but clearly distin- 
guishable curvature of the same sign and magnitude as predicted by our law 
(12. Hj) : hence we believe they have seen this law. 

The simulations were prompted by experimental DLCA studies carried 
out in the nineties by Earnshaw and Robinson, who among several other 
questions also tested O dj the validity of Aboav's law. Earnshaw and 
Robinson fitted their experimental m n versus 1/n data by a straight line and 
concluded [12] that DLCA accords well with the Aboav-Weaire law. We show 
one of their data sets in our figure El It exhibits the downward curvature 
that we expect. Our comment is again that whereas within the error bars 
a linear fit is also compatible with the data, only the curved theoretical law 
(12. 2|) provides an explanation. The experimental data seem only a small step 
away from being able to decisively rule out either one or the other. In view of 
the agreement observed by the authors between experimental and simulation 
data in many other respects [13], we are confident that DLCA experiments 
pushed to greater resolution will make the curvature visible. 

We finally ask what the physical reason is why RVF theory is applicable 
to DLCA. We believe the answer is as follows. One may obtain the Voronoi 
construction for randomly distributed seeds by starting a circular domain 
growth simultaneously from each seed and having it stop at the contact points 
of the circular frontiers. It is likely that in the experiments the depletion 
zones arising around each cluster play exactly the role of these growing circles. 
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Figure 5: Average sidedness m n of a neighbor of an n-sided cell. Full squares: 
DLCA simulation data by Fernandez- Toledano et al. (reference |34| . inset of fig- 
ure 5; the error bars were kindly provided to us by the authors). Open circles: 
experimental DLCA data by Earnshaw and Robinson (data for t = 60 min from 
|12j . figure 5). Thin straight lines: Aboav's linear law; the upper one is the fit to 
the simulation data and the lower one the fit to the experimental data proposed in 
reference [21] and in reference ^21 E]) respectively. Solid curve: first two terms of 
our theoretical large-n expansion for the Random Voronoi Froth, equation 12. 2|) . 
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Figure 6: Average number K\ of second neighbors of an n-sided cell. Full 
squares: DLCA simulation data by Fernandez- Toledano et al. (from [31], figure 
7; the error bars were kindly provided to us by the authors). Thin straight line: 
the linear fit to the data proposed in reference [21]. Solid curve: first term of our 
theoretical large-n expansion for the Random Voronoi Froth, equation H2.3|) with 
r] = 0.95. 
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3.3 Comparison to other exact models 

There are very few exact results for the two-cell correlation m n in other ge- 
ometrical systems. Those that we are aware of all concern model systems 
artificially constructed to be exactly solvable. An example is the anisotropic 
"laminated Poisson network" studied by Fortes [35] > for which Aboav's law 
appears to be satisfied exactly (for n > 4) with a = | and b = 14. Several 
results do exist, however, for topological planar cell models, i.e. models in 
which only the vertex connectivity is considered but where metric properties 
such as distances and areas are ignored. To this class belongs very interesting 
work by Le Caer [HE] and Le Caer and Delannay [22] • These authors start 
from regular planar ^-coordinated lattices with z > 3 and by a suitable algo- 
rithm construct from it a random topological froth that is three-coordinated 
just like the RVF and the vast majority of other cellular systems. For one of 
the simplest models of this kind (the case z = 4) Le Caer [3E] finds 

9 10 

m n = - + —, (n = 4,5,...,8), (3.2) 
2 n 

which is Aboav's law on a restricted n interval. Other topological models 
studied by these authors may be treated numerically exactly and, in the 
general case, show deviations from Aboav's law. 

Also in the class of topological models is work by Godreche et al. |3H] . 
who studied an ensemble of diagrams encountered in field theory. For these 
they found 
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m n = 7 + - + — — — , (3.3) 
n n[n + 1) 

which while curving upward approaches asymptotically the linear law nm n ~ 
In. 

It may well prove useful to classify models of cellular structure according 
to the decay of their m n by writing generically m n — ~ n~ a , where a is 
a positive exponent. Equations ()2.2|) and ()3.3|) then provide examples with 
a = | and a = 1, respectively. Beyond the RVF studied above it is unknown 
at present which other theoretical models and experimental systems have 
a — |, but for DLCA this certainly seems the likely value. We do not know 
if it is conceivable that a microscopic geometrical model of cellular structure, 
as opposed to the topological ones, could have a = 1. If so, its nm n would 
obey Aboav's linear law, either exactly for all n or asymptotically for large 
n. Obviously, knowing a system's a value is relevant for the analysis of its 
data. This value is very likely tied up with the cell formation mechanism, 
but in an as yet completely unknown way. 
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4 Beyond exact asymptotics: curve fitting 



Equation (|2.2|) is an exact expansion in powers of n 2 that we may write as 



1 



37T2 a 2 as , 
m n = 4+ — + ^ + 4 + (4.1) 

U2 n ri2 

The higher order coefficients a 2 , 03, ■ ■ ■ are well-defined and it is possible, in 
principle, to determine them successively. It would certainly be interesting 
to know a few more of them, but the corresponding calculations are far from 
straightforward. Moreover, it should be kept in mind that the asymptotic 
series (14. lj) is in all likelihood only asymptotic, which means (see e.g. |3Ti] . 
chapter 1) that even though the first few terms may give a good approxima- 
tion, when adding more and more terms for fixed n one eventually finds that 
the series does not converge. All we know is that as n grows, (J4.1j) approaches 
the exact result ever more closely. 

The remainder of this section is in a spirit different from the rest of our 
work. We include it because of the interest in numerical fits that exists in 
part of the cellular system community (see, e.g., reference jH]). In the same 
pragmatic way that led to Aboav's law we now ask if it is possible to fit 
the numerical m n data by a simple analytic expression that incorporates the 
first two terms of the asymptotic constraint ()2.2j) . The natural answer is to 
extend the two explicitly known terms of ()2.2|) by one or more extra terms 
with coefficients chosen to fit the data, and to truncate after that. With two 
extra terms one gets 



m r , 



4 + 3n^n~^ + An' 1 + Bn~%, (4.2) 



where A and B are on the same footing as Aboav's parameters a and b. With 
three extra terms one has 



11 1 _ 3 



m n = 4 + 37r2n~2 + An + Bn~z + Cn~. (4.3) 

Here A, B,C, . . . are not related in any simple way to a\, a 2 , 0,3, . . . of equation 
()4. lj) . A relation between the fit parameters arises if one imposes that the fit 
obey Weaire's sum rule [T7| 

(nm n ) = (n 2 ), (4.4) 

in which (X n ) = ^^=3 X n p n for any function X n of the sidedness n, and p n is 
the probability that an arbitrarily chosen cell have n sides. Thus, combining 
(I4.4j) with Aboav's law (jl.2|) leads to the well-known relation b = fi 2 + 6a 
where fi 2 = (n 2 ) — (n) 2 and one uses that (n) = 6. Imposing Weaire's sum 
rule makes good sense and we will therefore do the same. After substituting 
the assumed expressions ()4.2|) and ()4.3|) in ()4.4|) we find 

A + B{rT%) = 12 + /i 2 -37r5(n5). (4.5) 
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and 



A + B(rT*) + C(n- 1 } = 12 + /i 2 - 3vr5(n5), 



(4.6) 



respectively. The numerically most accurate data for p n and m n available 
today are those due to Brakke j2H]. They lead to ^ = 1-7807, (na) = 2.4344, 
(n~2) = 0.4160, and (n -1 ) = 0.1753, which are all needed in ()4.5|) and ([4.6)1 . 

We fitted both expressions ()4.2|) and (|4.3|) by minimizing the weighted 
maximum deviation from the data, i.e., letting denote Brakke's numerical 
values and a® their standard deviation, by searching for 



where the minimum is taken over the fit parameters under either constraint 
(14. 5|) or (|4.6j) . The results are shown in Table I. As witnessed by the value 
of e, the two-parameter fit still runs far outside of the standard deviations 
given by Brakke. The three-parameter fit yields strongly improved results. 
Comparison of the two- and three-parameter fit shows, in particular, that as 
more terms are added, there is no sign that the values of A, B, . . . converge. 
This also is a strong indication of the merely asymptotic character of the rT^ 
expansion. 

We emphasize that whereas sections El and El of our work rest on first 
principles, the procedure followed in the present section does not. It is hybrid 
in that it represents the construction of a best fit which has been made to 
include the first two terms of the exact large-n expansion. There is of course 
no reason to believe that m n can be represented as a sum of a finite number 
of powers of n, whether it be ()4.2j) . ()4.3j) . or Aboav's law (jl.2j) . 

5 Theory for non-RVF systems 

The RVF serves as a model of reference for general cellular systems much 
in the same way as an ideal gas does for interacting gases: some of the 
RVF result will apply to other systems and some will not. There are many 
factors that may potentially cause an experimental system to depart from 
RVF-like behavior. E.g., (i) in some cellular structures the 'seeds' represent 
actual particles or larger physical entities whose mutual interactions cannot 
be neglected; (ii) a general planar cell structure cannot be derived from a 
set of point centers by means of the Voronoi construction; and (iii) some 
systems, like soap froths, are not in equilibrium but rather in a coarsening 
state and need a dynamical theory in terms of cell transformation processes. 
We refer to Rivier |2Sj for an overview. 

We now ask the subtler question as to whether the cell-cell correlation 
m n in any of these non-RVF systems also deviates from Aboav's law by some 
slight curvature, and if it does, what the asymptotic behavior of its m n is. 



e = 



mm 



max \m n — m, 

3<n<15 




(4.7) 
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Two-parameter fit 


Three-parameter fit 


Simulation data 






A = 


3.816 




A = 


6.279 


(Brakke jggj) 






B = 


-7.163 




B = 


-18.327 














C = 


12.440 










50.46 




g = 


7.013 




n 




Tfl 

'"n 


777 — 777 




Tfl 

'"n 


777 — 777 


m B ± cr B 


3 


6 


9635 


-0 0481 
W . Ut:(j -L 


7 


01 83 


0067 


7 0116 ±0 0095 


4 


6 


7173 


0.0002 


6 


7152 


-0.0019 


6 7171 ±0 0028 


5 


6 


5005 


0.0083 


6 


4923 


0.0001 


6.4922 ±0.0017 


6 


6 


31943 


0.00442 


6 


31593 


0.00092 


6.31501 ±0.00014 


7 


6 


16816 


-0.00256 


6 


17114 


0.00042 


6.17072 ±0.00016 


8 


6 


04042 


-0.00974 


6 


04932 


-0.00084 


6.05016 ±0.00021 


9 


5 


93116 


-0.01622 


5 


94497 


-0.00242 


5.94739 ±0.00035 


10 


5 


83659 


-0.02200 


5 


85429 


-0.00431 


5.85860 ±0.00064 


11 


5 


7538 


-0.0254 


5 


7746 


-0.0047 


5.7793 ±0.0014 


12 


5 


6807 


-0.0297 


5 


7038 


-0.0066 


5.7104 ±0.0082 


13 


5 


6155 


-0.0364 


5 


6404 


-0.0115 


5.6519 ±0.0082 


14 


5 


557 


-0.0261 


5 


583 


0.000 


5.583 ±0.023 


15 


5 


504 


-0.0310 


5 


531 


-0.004 


5.535 ±0.070 



Table 1: Two-parameter fit ()4.2j) and three-parameter fit ()4.3j) to Brakke's 
simulation data |2H|- The latter are denoted here as m B and their statistical 
error as <r B . Both fits have been made to satisfy Weaire's sum rule (|4.4|) . 
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For certain microscopic models one may hope to be able to find the answer to 
this question perturbatively starting from the analysis of the present work. 

To show how such an approach might proceed we become more specula- 
tive. We discuss briefly and heuristically an example from class (ii) above, 
viz. the Voronoi tessellation associated with a gas of hard core particles of 
finite diameter d. For this system the preceding analysis of the large n limit 
remains valid as long as n < n*, where the crossover value n* is determined 
by the condition that the distance between adjacent first neighbors become 
comparable to the particle diameter. This gives 2nR c /n* ~ d whence, be- 
cause of R c = (n/47rA)2, it follows that n* ~ l/(Ac? 2 ), with the ~ sign 
indicating asymptotic proportionality. For n^n* the repulsion between the 
particles combined with the condition that they be locally aligned along a 
circle imposes that the radius of the central cell grows as R c ~ nd and hence 
that /s must saturate at a value f$ = c^X^d where Co is an unknown numeri- 
cal coefficient. The correlation m n can depend only on the two dimensionless 
variables n and Ac? 2 . It is reasonable to assume that in the scaling limit 
n — > oo, Xd 2 — > with nXd 2 = x fixed, m n may be expressed as 

m n ~ 4 + X^dM(nXd 2 ) (5.1) 

in which the scaling function Ai must satisfy 

37T 2 

M(oo) = c , M(x) ~ — j- for x — > 0. (5.2) 

£2 

Hence for hard core particles the limiting value of (|2.2jl is changed and 
()5.2j) does not tell us how this new limit is approached for large n. However, 
the n~2 decay law and the curvature effect that it entails survive in the 
crossover regime n < Ac? 2 , and this is the regime that is encountered first 
when the experimental or simulational precision increases. Future work will 
have to deal with this and other instances of deviations from RVF statistics. 



6 Conclusion 

Aboav's linear law (jl.2j) for the two-cell correlations m n has been known 
to fail for the Random Voronoi Froth. The exact calculation of this work 
has explained both qualitatively and quantitatively why it must fail: the 
exact asymptotic formula for m n exhibits an inverse square root decay with 
n. This large-n behavior represents a new paradigm in the field of planar 
cellular systems. 

In our discussion we have raised the question of whether similar violations 
occur also in other planar cellular systems, be they theoretical models or ex- 
perimental realizations. We found that departures from Aboav's law are 
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convincingly present in computer simulations of diffusion limited colloidal 
aggregation (DLCA) performed recently by Fernandez- Toledano et al. [34] . 
Since DLCA simulations have been found to be in all respects close to the 
corresponding experiments [14] . it is very likely that DLCA experiments vi- 
olate Aboav's linear law in the same way. Experiments by Earnshaw and 
Robinson [T21 EI] , m spite of error bars that do not unambiguously rule out 
Aboav's law, provide evidence for the curvature effect. One must expect that 
future experimental data (or perhaps even reanalysis of existing data), pro- 
vided they are of high enough precision and/or cover a large enough range 
in n, will reveal the presence of curvature in the m n versus 1/n relation also 
for certain cellular structures other than DLCA. 

The parameters a and b of Aboav's law are useful in that they provide a 
rough global classification of the behavior of a planar cellular system in the 
experimentally accessible range of measurement. No first-principle analysis 
of any geometrical cell model has, however, produced an interpretation and 
theoretical expressions for these parameters. Empirically, moreover, the pre- 
cise values of a and b depend on the range of the fit and on how it is performed 
[36J. Therefore in future high resolution simulations or experiments a fixa- 
tion on Aboav's law would be misguided and any curvature effect, whenever 
there is evidence for it, will be worthy of study. 
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ical period (CRCT) granted to the author by the French Ministry of Educa- 
tion in the academic year 2004-2005. The author thanks Dr. J.C. Fernandez- 
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K.A. Brakke for correspondence and Dr. N. Rivier for a discussion. 

A Appendix 

Equation ([2.2]) is the central result of this paper. Determining the coefficient 
of the second term in this equation involves a problem in statistics that we 
solve in this appendix. We have to prove the following theorem, illustrated 
by figure El 

Theorem. Let point particles ( 'seeds ') be randomly distributed with uniform 
density A in the upper half plane y > 0. Let Ti(x) be the piecewise parabolic 
curve that separates the upper half plane into a region of points closer to the 
x axis than to any of the seeds, and its complement. Consider the set of 
abscissae of the points where the parabolic segments join and where therefore 
Ti (x) has cusps. Then the density p of the abscissae on the x axis is equal 

to p = |Aa . 

Proof. The nontrivial part of this theorem is the coefficient §, since the 
seed density A may be scaled away; below we keep A only to have at all times 
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a check on the dimensionality of the quantities involved in the calculation. 
Let Pi = (xi,yi) be the position of the ith seed. The parabola 



= ! (A.i) 

Z V Hi 



separates the upper half plane into a region containing all points closer to 
the x axis than to Pi, and its complement. Each parabolic segment of Ti(x) 
lies on one of the fi(x). 

For the considerations that follow we must, as is always implied in sta- 
tistical mechanics, take a finite system and let eventually its size tend to 
infinity. We will start with a rectangular box [-L, L] x [0, L] whose volume 
we denote by V = 2L 2 [40J. Let the seeds inside this box be those of indices 
% = 1, 2, . . . , JV, where N is such that N/V = A. 

We choose a small Ax and define pAx as the probability that there is 
a cusp with abscissa in the interval [0, Ax]. This is the probability that 
there exist two parabolas, say y = fj(x) and y = fk{x), that have a point of 
intersection Pj k = (xjk, Vjk) with Xjk G [0, Ax] and such that Pjk is below all 
the parabolas y = fi(x) with £ ^ j,k. This may be expressed as 

pAx= f nlfiiT £ *^ a *) n w (a.2) 



where 



and 



l<i<N " u l<i<N l<j<k<N 1<£<N 



1 if < Xau < Ax 
Xjk (Ax) = { 3 (A.3) 

otherwise 



1 if f e (x jk )>y jk 

WVjk) = { _ . ( A - 4 ) 

otherwise. 



The sum in the integrand in (jA.2|) is on \N(N — 1) terms that all give 
identical results so that we may choose j,k = 1,2. The N — 2 integrals on 
(x 3 , y 3 ), ... , (xjy, Vn) ma y be carried out independently; we have 6^(2/12) = 1 
unless (x£,yt) is in the disk of radius y\ 2 centered at (xi2,yn)- Hence 

/L /*L 
dx £ dy e 9 e (y 12 ) = V -Tryl 2 (A.5) 
■L JO 



Using this we may rewrite (jA.2|) as 
pAx= — — ^ — / dxi / dx 2 / dyt dy 2 Xi 2 (Ax) [V - ny^] ~ 2 . 

J — L «/ — L J *j 

(A.6) 
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In the limit N, V, L -> oo with N/V = N/(2L 2 ) = A fixed this becomes 

/OO fOO pOO pOO 

dxi / dyi / dx 2 / dy 2 Xi2(Ax) exp(-A7n/ 2 2 ). (A. 7) 
oo J </-oo JO 



The exponential factor in the integrand of (jA.7|) has the interpretation of the 
probability that the disk of radius yi 2 around the point of intersection P 12 
contain no seeds. For i — 1, 2 we now transform from the pair of variables 
(xi, Hi) to the new pair (r*, s^) defined by 



r* = /*(0) = + 



2 7 ' 



Si = fm = -- i, (i = l,2). (A.8) 



The transformation has the Jacobian 



d(xi,yi) Ar-i 
In terms of these new variables of integration (jA.7|) becomes 



2 = 1,2). (A.9) 



pAx = 8A 2 / ds x / dri / ds 2 

JO J-oo JO 

dr2 (l + g 2p(i + g 2p Xw(Ax) exp(-A7n/ 2 2 ). (A.10) 

On the interval < x < Ax the two parabolas /j(x) may be represented by 
the expansion 

fi{x) = ri + Si x + 0(Ax 2 ), (2 = 1,2). (A.ll) 

The abscissa Xi 2 of their intersection point inside this interval, if there is one, 
is the solution of fifau) = j '2(^12)- Using (|A.11|) we find 

Xl2 = r -l^Il + 0(Ax 2 ). (A.12) 

S2 ~ Si 

The condition < x±2 < Ax imposed by X12 in ()A.10|) can be rewritten as 

n < r 2 < yi + (si - s 2 )Ax, (s!>s 2 ). (A.13) 

Because of symmetry we may restrict the s 2 integration in ()A.10|) to s 2 < s± if 
we compensate by an extra factor 2. Doing so, interchanging then the r 2 and 
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«2 integrals, using condition (|A.13|) . and still observing that yi 2 = r 2 + 0(Ax) 
we obtain from (|A.10I) 



pAx 



16A' 



dr 1 / dsi / ds 2 

J— oo J —oo 



J — oo 
fi + (si-S2)Ax 



dr 2 - 



r 1 r 2 



ri 



1 + ^) 2 (1 + ^ 



cxp (-Avrr 2 ), (A.14) 



valid up to corrections of order Ax 2 . Dividing by Ax and performing the r 2 
integral in the limit of Ax — > we find 



where 



p = I6\ 2 l r (h - I 2 ) 

/oo 
dr x r\ exp(— \-kt\] 
-oo 

= (47rAi)~ 1 



(A.15) 



and 



dsi / dss 



i + s ?) 2 (i + 4) 2 ' 



1,2). 



(A.16) 



(A.17) 



In the expression for I\ we interchange the s\ and s 2 integrals, whereas in 
the one for I 2 we interchange the names s\ and s 2 . The result is that 



h-h 



ds- 



*2 



A' 2 



Si dsi 



ds-, 



2(1 +s\: 



oo 




S2 


— oo/ 



ds 2 



2^3 



3tt 



(A.18) 



Upon combining ljA.15jl . (|A.16Jl . and (|A.18J) we find 

p = 16A 2 x (47rA§) _1 x (3tt/£ 



2 v ' 



(A.19) 



which is what we had to prove. □ 
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